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Abstract. We generalize Y. Nievergelt's inversion method for the Radon 
transform on lines in the 2-plane to the fc-plane Radon transform of continuous 
and L p functions on K n for all 1 < fc < n. 



1. Introduction 

Inversion formulas for Radon transforms of different kinds are of great impor- 
tance in mathematics and its applications; see, e.g., [^l4l[7l[8l[9l fl4l[l^[lTl[20l[24] . 
and references therein. Since many of them are pretty involved, especially for new- 
comers in the area, or applicable under essential restrictions, every "elementary" 
inversion method deserves special consideration. In 1986 Yves Nievergelt came up 
with intriguing paper [10] , entitled "Elementary inversion of Radon's transform" . 
His result can be stated as follows. 

Theorem 1.1. Let 

' l/(vra 2 ) if |i| < a, 



(1.1) G a (t) = 




VI - a 2 /* 2 



if \t\ > a; 



a > 0. Any continuous compactly supported function on the 2-plane can be recon- 
structed from the Radon transform over lines in this plane by the formula 



(1.2) 



1 f°° 

/(x,y) = lim — / / (Rf)(t — xcos a — y sina, a) G a (t) dt da, 
^ Jo J-oo 



where the double integral on the right-hand side equals the average of f over the 
disc of radius a centered at (x,y). 

Formulas (jl.ip and (|1.2p indeed look elementary. The following questions arise: 

1. What is the basic idea of the Nievergelt's method from the point of view of 
modern developments? 

2. Is this method applicable in the same elementary form to fc-plane Radon 
transforms on M ra for all 1 < k < n and arbitrary continuous or LP functions, 
satisfying minimal assumptions at infinity? 
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In the present article we answer these questions and indicate possible general- 
izations. 

Notation and main results. Let Q n ,h and G n ,k be the affine Grassmann manifold 
of all non-oriented fc-planes r in R™, and the ordinary Grassmann manifold of fc- 
dimensional subspaces C, of R n , respectively. Given £ G G n ,k, each vector x G R™ 
can be written as x = (x',x") — x' + x" where x' G C an d a; " G C^j C being 
the orthogonal complement to £ in R n . Each fc-plane r is parameterized by the 
pair (C,x") where ( G G n ,k and x" G The manifold will be endowed with 
the product measure dr = d^dx" , where dC is the S'0(n)-invariant measure on 
G ni fc of total mass 1, and dx" denotes the usual volume element on We write 
Co = Co(M. n ) for the space of continuous functions on R" vanishing at infinity; 
er„_i = 27r n / 2 /r(n/2) denotes the area of the unit sphere S n ~ 1 in R n . 

The fc-plane transform of a function / on R n is a function / on G n ,k defined by 

(1.3) /(r)= / f(x' + x")dx', T=(C,x")eG n 



k- 

c 

This expression is finite for all t G C/ nj fc if / is continuous and decays like 0(|a;| -A ) 
with A > k. Moreover [20l |22l [24] , if / G L P (R"), 1 < p < n/k, then /(r) is finite 
for almost all planes r G Gn,k- The above-mentioned bounds for A and p are best 
possible. 

Following [19] . we define the wavelet-like transform 



(1.4) {W* a y){x) = ± / ^(r) W (^-^)dr, a 



>0, 



where \x — r| denotes the Euclidean distance between the point x G R" and the 
fc-plane r. 

Theorem 1.2. [HI Th. 3.1] Let f e LP, 1 < p < n/k, and let ip(\ ■ |) be a 
radial function on R n , which has an integrable decreasing radial majorant. Lf w is 
a solution of the Abel type integral eguation 

r 

(1.5) cr 2 "" f 3 n - k - 1 w(s)(^-s 2 ) k / 2 - 1 d3^iP(r), c= a k-l °n-k-^ 



then 

(1-6) (W a */)0r) = y /(a - ay) V(|2/|) dtf, 

and therefore, 

(1.7) lira L(W a 7*)(a:) = A/(x), \= U{\x\)dx. 

a-,0 J 

R" 

T/ie Krmi m Jj. 7[ ) is understood in the L p -norm and in the almost everywhere sense. 
Lf f G Co H i p /or some 1 < p < n/k, then \1. 7| ) ZioZds uniformly on R ra . 

This theorem is a core of the convolution-backprojection method for the fc-plane 
Radon transform, and the most difficult task is to choose relatively simple functions 
w and ip satisfying (|1.5p . The crux is that the left-hand side of (|1.5j) has, in general, 
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a bad behavior when r — >oo. Hence, to achieve integrability of if), the solution w 
must be sign- changing. 

Our first observation is that the essence of Y. Nievergelt's Theorem 11.11 can be 
presented in the language of Theorem 11.21 as follows. 

Theorem 1.3. Let k = 1, n = 2, 

(1.8) ^(r) = 

Then il.5\) has a solution 



1 if < r < 1, 
if r > 1. 



(1.9) w{r) = 



1 if < r < 1, 



1 ; if r > 1. 



smc/i i/iai /or every compactly supported continuous function f on R 2 , 

(1.10) (WC /)(*)= / f(x-ay)dy, 

\v\<i 

where W* is the wavelet-like transform generated by w. 

For the convenience of presentation, we will keep to the following convention. 

Definition 1.4. The convolution-backprojection algorithm in Theorem 1 1.21 will be 
called the Nievergelt 's method if ip is chosen according to (|1.8p . 

Of course, Theorem 11.21 deals with essentially more general classes of functions 
than Theorem 11.11 however, the main focus of our article is different: we want to 
find auxiliary functions tp arid w, having possibly simple analytic expression. 

Theorem A. 

(i) The Nievergelt's method is applicable to the X-ray transform (the case k — 1 ) in 
any dimension. Namely, if 'ip is chosen according to A1.8\) , then U.5\) has a solution 

(I if < r < 1, 



(1.11) w(r) = < 



r 3-n 1 



2v^Fr 



v n/2-l ( r 2_ w )-3/2 dw if r> l ) 



and inversion formula JJ.7| ) holds. 

(ii) Ifk>\, then Nievergelt's method is inapplicable. 

An integral in (|1.11[) can be expressed through the hypergeometric function and 
explicitly evaluated in some particular cases; see Section 2.2. For instance, if n = 2, 
then p. lip is the Nievergelt's function (|1.9p . 

To include all 1 < k < n, we modify the Nievergelt's method by choosing ip(r) 
in a different way as follows. 

Theorem B. Let 

if < r < 1, 

(1-12) ^)=<| 

r n + 2< + l % f " > V > 
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I > 0. Then the corresponding function tp(\ ■ |) has a decreasing radial majorant in 
L (R n ) and hl.5\) has the following solution: 



(i) In the case k = 2£; £ = 1,2,.. 



(1.13) w(r) = 



2+21-n 



c e r 



(ii) In the case k = 21 + 1; £ = 0, 1, 2, 


(1.14) w(r) 



(--) 


'(r 2 -l) r 


\2r dr J 


r 2«+3 


r(n/2- 






i/ < r < 1, 
if r > 1, 



r(n/2) 



3+2f-n 



2r rfr 



^ r 2_ 1 ^£+l/2 



r 2£+2 



«/ 0<r<l, 



if r > 1, 



T((n-l)/2-£)£\ 



T(n/2)T(£ + 3/2) ' 
In both cases the inversion result in Theorem \1.2\ is valid. 

Theorems A and B are proved in Sections 2 and 3, respectively. 

Possible generalizations. 

1°. The convolution-backprojection method is well-developed in the general context 
of totally geodesic Radon transforms on spaces of constant curvature. Apart of W 1 , 
the latter include the n-dimensional unit sphere S n and the hyperbolic space H n ; 
see [TJ [TH1 HH] . As above, the key role in this theory belongs to a certain Abel type 
integral equation and the relevant sign-changing solution w. Moreover, passage 
to the limit in (|1.7j) as a — * 0, can be replaced by integration in a from to 
oo against the dilation-invariant measure da/a. This leads to inversion formulas, 
which resemble the classical Calderon's identity for continuous wavelet transforms 
[3]. The corresponding wavelet function is determined as a solution of a similar 
Abel type integral equation; see [TJ [THl [19] for details. In all these cases analogues 
of Theorems A and B can be obtained. We leave this exercise to the interested 
reader. 

2°. Unlike the classical fc-plane transforms on K n , the corresponding transforms on 
matrix spaces [SJ [TTJ [T^l [T3J are much less investigated. To the best of our knowl- 
edge, no pointwise inversion formulas (i.e., those, that do not contain operations in 
the sense of distributions) are available for these transforms if the latter are applied 
to arbitrary continuous or IP functions (p ^ 2). One of the reasons of our interest 
in Nievergelt's idea is that it might be applicable to the matrix case. Moreover, 
as in 1°, it may pave the way to implementation of wavelet-like transforms in the 
corresponding reconstruction formulas. We plan to study these questions in our 
forthcoming publication. 
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2. The case fc = 1 

2.1. Proof of Theorem A. We will be dealing with Riemann-Liouville fractional 
integrals 

1 f u 

(2.1) (C + g)(u) = — (u-vr^g^dv, a>0. 

J a 

Changing variables, we transform the basic integral equation 

r 

(2.2) cr 2 ~ n f ^-^wis^r 2 - s 2 ) fe /2-i ds = ^(r), c ^ a k-KT n -k-i 



C~n-1 


(cf. (HHl)) to the form 

'0- 



(2.3) I^ 2 w = 



where 

9 n/2-1 

(2.4) w(u) = u^/^wi^l), m = 

Suppose that ip is defined by (|1.8p . If < ti < 1, then, by homogeneity, 

(2.5) w(u)-CU , c - cr(fc/2)r((n _ fc)/2) 

c being the constant from (|2.2j) . Hence, for u > 1, we necessarily have 

(2.6) (I^^H = J\u-v^ 2 -h^ 2 -'dv. 

If fc > 2, this equation has no solution w G Lj oc (l, oo), because, otherwise, we get 
lim (l.h.s) = 0. lim (r.h.s) = const ^ 0. 

This proves the second statement in Theorem A. 

Consider the case fc = 1. If < U < 1, then m}(m) = i/^™ -1 )/ 2-1 . If u > 1, then, 
setting u>o(u) = u( n '~ 1 )/ 2_1 , from (|2.6[) we get 

(^{ a [t«-t«o])(u) = p / (u-w)- 1 / 2 !)^- 1 )/ 2 - 1 ^ 



_ /2 _ t r((n-l)/2) 

T(n/2) 



This gives 



w(u) = ^(u) - ^= -1 /""(« - i;)- 1 / 2 !;"/ 2 - 1 d V 



rw du 



or (set f" = - Jo) 



?i ( M ) = _^ / (u-t,)- 3 / 2 ^/ 2 - 1 ^. 
1 ^ 2^ io 



Thus, if Iq+w = rp, then, necessarily, 



(2.7) u)(u) 



z (n-l)/2-l ifO<M<l, 
el 

(u-tj)- 3 / 2 ^"/ 2 -^^ ifu>l. 



2v^ Jo 
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One can readily see that function (|2.7|) is locally integrable on (0, oo). Let us 

1/2 ~ 

prove that it satisfies (I ^_ w)(u) = ^(u) for all u > 0. It suffices to show that 
(Io+w)(u) = when u > 1. We have Iq+w = I\ — 1%, where 



Jo 



h = — / (u-v)-^ 2 dv (v-sy^s^-Us. 

Both integrals can be expressed in terms of hypergeometric functions. For I\ , owing 
to 3.197 (3) and 9.131 (1) from [BJ, we obtain 

V 2 ' J \2 '2 ' 2 ' u 
r((n-l)/2) 1/2 / l.n+1, 1 

(2 - 8) - V^r((n + l)/2) {U - 1] F { 1 '2'^'— U/ 

For I 2 , changing the order of integration and using [BJ 3.238 (3)], we have 



T c„ (u - l) 1 / 2 f 1 s"/ 2 - 1 (1 - aJ-Va 

i 2 = / 



ds. 

s 

By [BJ 3.228 (3)] this expression coincides with (|2.8| . 

To complete the proof, we recall that w(r) = r 3 ~ n w(r 2 ), which gives 

if < r < 1, 

3-n r 1 



w(r) 

i *- 

This coincides with 



( r 2 _ ^-3/2^/2-1^ if r>1 . 



2.2. Examples. Let us give some examples of functions w(r) defined by (jl.lip in 
the case r > 1. By [BJ 3.197(3)], 

l j 2^Fr(n/2 + l) V2 2' 2 ' r 2 

Keeping in mind that -F(a, 6; c; z) = F(b, a; c; z), and using formulas 156, 203, and 
211 from 15, 7.3.2], we obtain: 

For n = 2: 



For n = 3: 



1 ^3 3 5 1 \ 2/ .1 1 

= 1 ^ o? o! — J=-(arcsm 



For n = 4: 



2r 3 0Fr(5/2) V2'2'2'r 2 y tt V r V^Tf 

... _ , ... . 2r 2 - 1 
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3. The general case 
As in Section |2~T1 our main concern is integral equation (|2.2j) . which is equivalent 



to 

(3.1) 
(3.2) 



T k/2 ~ 7 

I + w = ip, 

5(«) = vl n - k)/ *- 1 w(,Vu), ^(«) = 2 ""f/J> ^(V")- 



c r(fc/2) 

We want to find relatively simple functions w and ?/>, which are admissible in the 
basic Theorem 11.21 and such that the corresponding functions w and ip obey (|3.1[) . 
It is convenient to consider the cases of k even and A; odd separately. 

3.1. The case of k even. Let k — 2£; £ = 1,2, We choose 

if < r < 1, 



(3.3) 



ip(r) = < 



(r 2 - 1)* 



if r > 1. 



The corresponding function x — > VKI 3 -!) obviously has a decreasing radial majorant 
in i^R"). By (f3~2j) . equation (f5TT|l is equivalent to 



(3.4) (/*+«>)(«) 
The £th derivative 



c e (u - 1)* 



^+3/2 



r(n/2-^) 
r(n/2) ' 



u > 1. 



(3.5) w(u) 



Cj (u - I) 1 
u e+3/2 



u-lf] 



is integrable on (l,oo). The fact that (|3 . 5[) satisfies (|3.4[) can be easily checked 
using integration by parts. Thus, the pair of functions wir) and ip(r), defined by 



(3.6) w(r) 



ct r 



2+21-n 



(--)' 


"(r 2 -l) r 


\2r dr J 


r 2t+Z 



if < r < 1, 
if r > 1, 



and (|3.3p . falls into the scope of Theorem 11.21 and the "even part" of Theorem B is 
proved. 

3.2. The case of k odd. Let k = 2£ + 1; £ = 0, 1, 2, . . . . We define ?/>(r) by lpT3]) . 
as above. Then, instead of (13.41). we have 



(3.7) (lit 1/2 w)(u) 



u e+3/2 



T((n -l)/2-£) 
r(n/2) • 



71 > 1. 



This gives = 5 < - £+1 - > (u), where 



9{u) = ~j= 



(n-u)- 1 / 2 ^-!^ _ ( u _ 1 )^+i/2 
— dv = ci - 



ci = 



v e+3/2 

T((n- l)/2-l)l! 
r(n/2)r(£ + 3/2) ' 



s 
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(use [6j 3.238(3)]). Let us show that w(u) = g^ +1 \u) satisfies (|3.7|) . Integrating 
by parts, we have lft 1 ^ 2 'W = l\+g', where 



^ dg(u) 
9 (u) = —, — = c t 
du 



(£ + 1/2) (u - if- 1 ' 2 (£ + 1) (u - l) l+1 ' 2 



,e+i 



,1+2 



This gives I 1+ 1 w = q (Zi — I 2 ) where (use [51 3.238(3)] again 
/i = 



£+1/2 f (u- v)- 1 ' 2 ^- l)^ 1 / 2 , r(£ + 3/2) (u-lV 

dv = 



f.+i 



/+1/2 



h = 



C+l f (w- U )- 1 /2( l ,_ 1 )«+i/2 r(£ + 3/2) (u- 1) £+1 

aw = 



,4+2 



,£+3/2 



Hence, 



and therefore, 



h-h 



r(£ + 3/2) (u-iy 



i\ u e + 3 / 2 



Ui+ w)(m) 



5,r(£ + 3/2) ( U -l) £ _ , ( U -l) £ 



£! u«+3/2 
as desired. Thus, functions w(r) and tp(r), defined by 





£+3/2 



V 



(3.8) w{r) 



a r 



3+2£-n 



}_d_\ +1 \{r 2 - l) e+1/2 
2r dr 



r 2£+2 



if < r < 1, 



if r > 1, 



and (|3.3p . obey Theorem 1 1.21 This completes the proof of Theorem B. 
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